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Abstract 

A  method  for  remote  optical  measurements  of  magnetic  field  variations  above  the  surface  of 
seawater  is  analyzed.  This  magnetometry  mechanism  is  based  on  the  fact  that  polarized  laser 
light,  in  the  presence  of  the  Earth’s  magnetic  field,  will  exhibit  a  polarization  rotation  when 
reflected  off  the  surface  of  the  water  and  off  an  underwater  object.  The  two  mechanisms 
responsible  for  the  polarization  rotation  are  the  Surface  Magneto-Optical  Kerr  Effect  (SMOKE) 
and  the  Faraday  effect.  In  both  mechanisms,  the  degree  of  polarization  rotation  is  proportional  to 
Earth’s  local  magnetic  field.  Variations  in  the  Earth’s  magnetic  field  due  to  an  underwater  object 
will  result  in  variations  in  the  polarization  rotation  of  the  laser  light  reflected  off  the  water’s 
surface  (SMOKE)  and  off  the  underwater  object  itself  (Faraday  effect).  An  analytical  expression 
is  obtained  for  the  polarization-rotated  field  when  the  incident  plane  wave  is  at  arbitrary  angle 
and  polarization  with  respect  to  the  water’s  surface. 
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I.  Introduction 

Optical  magnetometry  can  be  a  highly  sensitive  method  for  measuring  small 
variations  in  magnetic  fields  [1-3].  The  development  of  a  remote  optical  magnetometry  system 
would  have  important  applications  for  the  detection  of  underwater  and  underground  objects  that 
perturb  the  local  ambient  magnetic  field.  For  a  number  of  magnetic  anomaly  detection  (MAD) 
applications,  juG  magnetic  field  variations  must  be  detected  at  standoff  distances  greater  than 

lO^mfrom  the  sensor  [4-5].  The  remote  atmospheric  optical  magnetometry  mechanism 
considered  here  is  based  on  polarization  changes  in  reflected  laser  light  from  the  surface  of  sea 
water. 

The  dielectric  properties  of  water  in  the  presence  of  a  magnetic  field  will  rotate  the 
polarization  of  an  optical  field.  The  polarization  rotation  occurs  via  two  distinct  mechanisms: 
the  Surface  Magneto-Optical  Kerr  Effect  (SMOKE)  and  the  Faraday  rotation  effect.  SMOKE  is 
purely  a  surface  phenomenon  in  which  the  degree  of  polarization  rotation  of  the  reflected  light  is 
proportional  to  the  magnetic  field  but  not  to  the  propagation  distance  in  the  water  [6].  The 
SMOKE  mechanism  is  used  in  material  science  research  in  such  devices  as  the  Kerr  microscope. 
The  Faraday  effect  is  due  to  the  difference  in  propagation  speed  of  right-hand  and  left-hand 
polarized  light,  and  the  resulting  degree  of  polarization  rotation  is  proportional  to  both,  the 
magnetic  field  and  the  propagation  distance. 

In  this  paper,  we  consider  whether  these  effects  can  be  used  for  the  remote  detection 
of  underwater  objects.  An  underwater  object  can  produce  appreciable  magnetic  field 
perturbations  above  the  water’s  surface.  The  proposed  magnetometry  mechanism  is  shown 
schematically  in  Fig.  1.  Here,  a  linearly-polarized  laser  beam  is  propagated  to  and  reflects  off 
the  water’s  surface  near  the  detection  site.  The  degree  of  polarization  rotation  in  the  reflected 
field  is  proportional  to  the  local  magnetic  field.  The  magnitude  of  the  polarization-rotated 
component  of  the  reflected  field  is  measured  at  various  locations  near  the  inspection  site.  An 
irregularity  in  the  measurement  at  one  location  would  correspond  to  a  deviation  in  the  local 
magnetic  field.  The  local  field  deviation  could  indicate  the  presence  of  an  underwater  object. 
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Figure  1 :  Polarized  laser  light  is  propagated  to  the  surface  of  the  water  and  exhibits  polarization 
rotation  upon  reflection  due  to  Faraday  rotation  and  SMOKE.  The  purpose  of  this  paper  is  to 
estimate  the  polarization  rotation  angle  Ad  cc  .  If  A6>  is  large  enough  to  be  precisely 
measured,  the  mechanism  shown  in  this  figure  enables  the  measurement  of  fluctuations  in  the 
Earth’s  magnetic  field  ,  which  in  turn  may  indicate  the  presence  of  an  underwater  object. 
The  solid  arrows  in  the  figure  denote  wave  vectors. 


II.  Model 

To  analyze  the  Surface  Magneto-Optical  Kerr  and  the  Faraday  rotation  effects,  we 
first  obtain  the  dielectric  properties  of  magnetized,  pure  water.  The  polarization  field  associated 
with  magnetized  water  at  the  laser  frequency  is  obtained  using  an  extended  form  of  the  Lorentz 
model  [7].  The  polarization  field  is  used  to  derive  a  wave  equation  for  the  reflected  laser  field. 
By  appropriately  matching  boundary  conditions  at  the  various  interfaces,  we  obtain  the  reflected 
field  which  has  its  polarization  modified  by  both,  the  SMOKE  and  Faraday  effect.  The 
polarization-rotated  component  of  the  field  is  proportional  to  the  local  magnetic  field.  In 
applying  the  boundary  conditions,  the  laser  propagation  direction  and  polarization  as  well  as  the 
Earth’s  magnetic  field  are  taken  to  be  in  an  arbitrary  direction. 
i)  Polarization  Field  of  Magnetized  Water 

In  the  Lorentz  model,  the  dielectric  properties  of  the  molecule  are  represented  by  the 
motion  of  the  forced  displacement  of  the  electron  distribution  from  its  equilibrium  position.  The 
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displacement  of  the  electron  distribution  5r{Y,t)  is  modeled  by  the  following  driven  harmonic 
oscillator  equation, 


5y 


E+^xB  Ic  , 
dt 


(1) 


where  Qg  is  the  electron  binding  frequency,  y  is  the  damping  coefficient,  q  is  the  electronic 
charge,  m  is  the  electronic  mass,  andB^  is  the  local  magnetic  field,  i.e..  Earth’s  magnetic  field 

~  0.5  G.  The  laser  electric  field  is  E(r,r)  =  Re(^E(r)exp(-/ft>r)^  where®  is  the  frequency. 

The  magnetic  field  associated  with  the  laser  is  neglected  on  the  right  hand  side  of  Eq.  (1)  since  it 
produces  forces  at  frequencies  of  2®  and  0. 

The  polarization  field  associated  with  the  magnetized  water  molecules  is  given  by 

P  (r)  =  qNdf ,  where  df{r)  represents  the  electronic  displacement  from  equilibrium  of  the 

water  molecule,  N  is  the  effective  bound  electron  density  of  the  water  molecules,  and  phasor 
notation  is  used  throughout  to  represent  the  various  field  quantities,  i.e., 

Q{r,t)  =  Re(0(r)  exp(-/®r)^.  The  displacement  5f{r)  is  given  by 

-  iYO)^Sf(r)  =  q{^E(r)  -  io)Sf(r)xB^  / c^/m.  (2) 

The  polarization  field,  correct  to  first  order  in  the  magnetic  field,  is 

2 

P(r)  =  qNSf^  y  (E(r)  +  io,n,xE(r)/ D{m)),  (3) 

47rDy(Oj '  ^ 

where  Di^co)  =  -  oP  -  iyco ,  =  {ATzq^N /  ,  and  £!„  =  qB^^  / me  ~\Qi’ s~^ . 

The  laser  frequency  is  chosen  to  be  in  the  blue-green  range  where  the  penetration  depth  is 
maximum.  In  the  following  examples,  the  laser  frequency  is  taken  to  be  co  =  3.8xl0'^rad/s 
(®  =  2.48eV /  /z)  where  pure  water’s  absorption  coefficient  is  minimum,  i.e.,  a  =  5x10“^ cm  ' ; 
the  corresponding  penetration  depth  (e- folding  length)  is  /z^  »  20  m  .  In  the  blue-green  regime, 
the  real  part  of  the  refractive  index  is  Re[zZo] « 1.34 .  The  lowest  energy  electron-binding 
frequency  in  water  is  Qg  «  23.5eV /  h ,  »  21eV  /  fi ,  and  the  damping  coefficient  is 

y  «3x10‘^5'“'. 
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In  general,  there  are  many  resonances  above  which  contribute  to  the  refractive  index  of 
water.  We  neglect  these  binding  frequencies  since  they  are  at  higher  energies.  In  addition,  the 
oscillator  strength  parameter  associated  with  the  resonance  is  taken  to  be  unity  [8]. 

Nonlinear  effects  in  water  can  be  neglected  since  the  laser  intensity  is  well  below  the  level  where 
these  become  important. 
ii)  Wave  Equation 

In  the  absence  of  free  charges  and  currents.  Maxwell’s  equations  for  the  fields  in 

/V  /V  /V  /V  /V  /V 

phasor  notation  are  given  by  VxH  =  -/(ty/c)D,  V-D  =  0 ,  VxE  =  i(ft)/c)B  andV-B  =  0. 
These  equations  can  be  combined  to  give 

V"E  +  (^u"/c")D  =  V(V-E).  (4) 

To  obtain  a  wave  equation  for  E  ,  we  note  that  the  electric  flux  density  D  =  (E  +  4;tP)  is  given 

/\  /  \ 

by  D  =  +  i{s^-\){col D{co))Ei.^  X  )E,  where  the  electric  polarization  field  is  given  in  Eq.  (3) 

and  where  the  dielectric  constant  of  water,  in  the  absence  of  the  Earth’s  magnetic  field,  is 
{of)  =  1  +  /  D{of) .  Since  V  •  D  =  0 ,  we  find  that 

=  ”  V.(n,xE)  =  3^°~‘>  (n.  (VxE)),  (5) 

D{af)  D{af) 

where  the  final  form  for  V  •  E  is  valid  since  is  spatially  uniform.  Using  Eqs.  (4)  and  (5),  the 
wave  equation  takes  the  form 


^  CO  CO  CO  (s  — D  CO 

72^  .  ^  -(Q^xE) +  ’  V(Q„-(VxE)), 


V^E  +  ^„^E  =  -/^(^,-l) 


"  '  D{of) "  "  "  D{of) 

The  right-hand  side  of  Eq.  (6)  can  be  simplified  by  noting  that 

V(Qp  •  (V  X  E))  =  X  (V  X  V  X  E))  -i-  •  V)(V  x  E)  and,  to  lowest  order  in  , 

V  X  V  X  E  =  {o)^  /  c^)£’^E  .  The  final  form  for  the  wave  equation,  correct  to  first  order  in  ,  is 

,{£,-l)  0) 


(6) 


V  *^7 


E  =  /- 


-(a„-v)(VxE). 


(7) 


D{of) 

Since  the  right-hand  side  of  Eq.  (7)  is  small,  we  can  write  the  total  field  as  E(r)  =  E^  +  SE, 
where  E^  here  denotes  the  laser  field  in  the  absence  of  the  perturbation  (i.e.  the  Earth’s  magnetic 
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field)  and  SE  is  the  small  polarization-modified  component  of  the  laser  field  due  to  the  Earth’s 
magnetic  field  and  driven  by  the  laser  field  .  Separate  wave  equations  can  be  written  for 
these  fields: 


r 


,2  A 


V  ^  J 


E  =0, 


V  ^  J 


^  ^  (a„-v)(VxEj. 


Dioj) 

In  the  special  case  where  the  laser  field  E^  is  transverse  to  the  magnetic  field  ,  one  can 
recover  from  Eq.  (8b)  the  usual  Faraday  rotation  angle,  i.e. 

Qq  L/(2c  Qg)  □  (10^G”'m”'rad)Lfi^ ,  where  L  is  the  propagation  distance  [9]. 

These  wave  equations  are  numerically  solved  for  arbitrary  laser  field  propagation  direction, 
polarization  and  magnetic  field  direction. 


(8a) 

(8b) 


III.  Analysis  of  Polarization  Rotation 

The  proposed  magnetometry  method  is  feasible  only  if  the  laser  beam’s  polarization 
rotates  a  substantial  amount  upon  reflection  off  the  water’s  surface.  The  purpose  of  the  analysis 
is  thus  to  calculate  the  amount  of  polarization  rotation  A^=|  JE’l/lEgjI,  where  (JE  is  the 

polarization-rotated  component  of  the  reflected  laser  field  and  is  the  un-rotated  component. 

Alternatively,  this  is  written  as  A6'=(l  /R)|<5’£'|/|£'o|  ,  where  is  the  incident  laser  field  and  the 
reflection  coefficient  is  R  =  (n  -  l)/(n  -l- 1)  «  0. 1 5 . 

We  first  consider  a  simplified  configuration  as  shown  in  Fig.  2  to  obtain  estimates 
for  the  amount  of  polarization  rotation  A  <9 .  Here  a  linearly-polarized,  monochromatic  laser 
beam  is  propagated  parallel  to  Earth’s  magnetic  field  which  is  taken  to  be  perpendicular  to  the 
surface  of  the  water.  An  optically  reflective  underwater  object  at  depth  h  below  the  water’s 
surface  is  represented  by  a  perfectly  conducting  plate  parallel  to  the  surface. 

From  this  simplified  configuration  will  emerge  two  distinct  contributions  to  the 
polarization  rotation:  the  Faraday  effect  and  SMOKE.  The  Faraday  effect  yields  the  ratio 
|<9£'/£'o|~12.7rnQ ,  whereas  SMOKE  yields  3x10“'^ .  The  Faraday  effect  is 

larger  by  the  ratio  h/  AU  1  . 
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The  SMOKE  contribution  to  the  polarization  rotation  is  obtained  for  an  arbitrary 
configuration. 
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Figure  2:  A  simple,  geometrical  configuration  of  reflected  laser  light  undergoing  polarization 
rotation  due  to  Faraday  rotation  and  SMOKE.  The  purpose  of  this  configuration  is  to  obtain 
estimates  for  the  polarization  rotation  angle  /S.0 . 


The  solid  arrows  in  Fig.  2  are  the  wave  vectors  while  the  dashed  arrows  denote  the  electric  fields 
corresponding  to  the  incident  wave  ,  the  transmitted  wave  E,. ,  the  object-reflected  wave  E^2  > 

and  the  surface-reflected  wave  E^j  which  is  modified  by  a  small  polarization-rotated  component 
S'E .  The  electric  field  can  be  written  as, 

E{z)  =  E^Qx^{ikz)e^+Ej^^Qx^i-ikz)e^+dEQx^i-ikz)ey  ,  (z<0)  (9a) 

E (z)  =  E^  ex^{ik^z)e^  +  E~  ex^{ik_z)e_  + 

£■^2  exp(-il^z)e^  -t  £^2  exp(-i£_z)e_ ,  (0  <  z  <  /?) 

where  ={e^±iey)  /  ^/2  ,  and  are  unit  vectors  in  the  x  and  y  direction  respectively, 

A:  =  ro/  c  is  the  wavenumber  in  air,  =ck+  /  co=^s^±s-^  is  the  complex  index  of  refraction  in 

magnetized  water,  and  =  {e^  -  1)  co€ll  D{co) .  Application  of  the  electromagnetic  boundary 

conditions  yields  for  the  configuration  shown  in  Fig.  2  the  following  ratio  for  the  polarization- 
rotated  field. 
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SE 


n.  cos 


(Kh) 


n  cos 


{kh) 


cos{k^h)—i?,m[k^h)  n_  cos(^_/?)-/sin(A:_/?)^ 

In  the  limit  that  /^  ^  co ,  the  refleeted  Faraday  rotation  effect  vanishes  and  the  polarization 
rotation  is  due  solely  to  SMOKE.  In  this  limit  Eq.  (10)  reduces  to 

SE 


(10) 


-n_ 


^0  (l  +  «+)(l  +  «-) 

where  ^  -l)(co  Qg  /  Q^)  /  (n^  +ng)  □  3  x  10“'^ .  The  ratio  of  the  fields  given  in  Eq.  (10)  is 
shown  in  Fig.  3  as  a  function  of  the  depth  of  the  eonducting  plate  h  . 


(11) 


h  (meters) 


Figure  3:  Ratio  of  the  polarization-rotated  reflected  field  to  the  ineident  field  as  a  fiinetion  of 
object  depth  h  .  The  polarization  rotation  results  from  a  eombination  of  Faraday  rotation  and 
SMOKE;  however,  the  contribution  due  to  SMOKE  is  small. 


We  now  consider  the  general  SMOKE  contribution  to  the  polarization  rotation  by 
analyzing  the  proeess  for  arbitrary  field  direetions.  The  polarization-rotated  field  due  to  SMOKE 
as  a  function  of  the  field  components  of  an  ineident  plane  wave,  for  arbitrary  polarizations  and 
angle  of  ineidence  as  well  as  arbitrary  orientation  of  the  magnetic  field,  is  obtained  in  Appendix 
A  and  found  to  be, 

|^£|  □  K{e)|n,|(£;i-”  +2{ElE]^+E':-El)c<xf  e+ElE';  cos-  0^ ,  (13) 
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where  K(^)  «  (Qq  co !  g{0)  =  2.5x10“'^  g{6)  □  2.5x10“'^  and 

rig  «  (cos^g  -  6>sin^  6g  sin^  (pgll.K)  ~\  .  In  Eq.  (12),  6  is  the  angle  of  the  incident  wave  vector 
with  respect  to  the  water’s  surface  normal  vector.  The  angles  6g  and  are,  respectively,  the 

polar  and  azimuthal  angle  of  Earth’s  magnetic  field  with  respect  to  the  incident  wave  vector. 
Appendix  A. 


g(^) 


Figure  4:  The  function  g{9)  plotted  above  represents  the  angular  dependence  of  the 
polarization-rotated  field  on  the  incident  wave  vector’s  orientation  for  a  specific  geometry 
{0g=rl>g=Q  ,  E^=Q). 

We  find  that  no  choice  of  polarization  or  angle  of  incidence  will  significantly 
increase  the  SMOKE  contribution  to  the  polarization  rotation:  see  Fig.  4. 

IV.  Discussion 

A  mechanism  for  the  remote  optical  detection  of  a  shallow,  i.e.,  ~  1  -  10  m, 
underwater  conducting  object  has  been  presented.  In  this  mechanism,  a  polarized  blue-green 
laser  beam  is  reflected  off  the  water’s  surface  and  the  reflected,  polarization-rotated  electric  field 
recorded.  A  variation  in  the  measurements  at  two  nearby  locations  may  indicate  the  presence  of 
an  underwater  conducting  object.  The  polarization  rotation  is  mainly  due  to  the  Faraday  effect; 
SMOKE  has  a  negligible  contribution. 

As  an  alternative  to  the  mechanism  proposed,  one  can  detect  the  underwater  object 
directly  by  measuring  the  un-rotated  component  of  the  reflected  laser  field.  This  reflected 
component  can  be  significantly  larger  in  amplitude  than  the  polarization-rotated  component.  By 
time  gating  the  laser  pulses,  the  object’s  depth  can  be  determined. 
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Appendix  A 

A  perturbative  solution  is  described  here  for  the  polarization-rotated  reflected  field 
associated  with  SMOKE.  The  treatment  allows  for  the  incident  plane  wave  to  have  arbitrary 
angle  and  polarization  with  respect  to  the  water’s  surface  and  allows  for  arbitrary  magnetic  field 
orientation.  The  relevant  coordinate  systems  are  labeled  in  Fig.  Al. 


Figure  Al :  Coordinate  systems  used  in  the  calculation  of  the  polarization-rotated  reflected  field 
associated  with  SMOKE  given  an  arbitrary  geometric  configuration. 


The  arrows  correspond  to  the  transmitted  wave  vector  in  the  absence  of  a  local  magnetic  field 
and  the  incident  wave  vector;  the  boundary  conditions  ensure  that  they  are  coplanar  and  that 
Uq  sin  a  =  i^co/c)  sin  6  ,  where  is  the  magnitude  of  the  transmitted  wave  vector.  The 
coordinate  systems  are  defined  such  that  x  =  x'  =  x”  coincide. 

The  fields  E(r)  ,  D(r)  ,  B(r)=-i(c  /  co)VxE(r)  are  expanded  to  first  order  in  the 
small  parameter  A=a)ClQ  / D(ci>)  □  10“"  where,  as  before,  D[ol>)  =  -  0/  -  iyco  .  Since 

V-.D(r)=0,  to  each  order  in  A  ,  Eq.  (4)  yields  the  wave  equation  V^E^  +SQ{oy^  I c^)Eq  =0, 

where  the  field  subscript  denotes  zeroth  order  in  A  .  This  wave  equation  is  used  to  calculate  the 
zeroth  order  field  amplitudes. 

For  z'  >  0 ,  Eq.  (7)  implies  that,  for  the  transmitted  field Ej,(r)=^y  ^(r)  -1-  dEjir) , 

+0/  £q  / c^^SEj,(r)=i(^£Q-l)co(Q.^^  x  ^(r)^ (^0  D(ct>))  (13) 

The  field  is  rewritten  as  (5'£y(r)=r(r)exp(/A:oz")=r(r)epexp(iA:oz")  in  order  to  employ  the 
paraxial  approximation  -^-k^  +2ikQd / dz” ,  which  states  that  the  field  varies  a  negligible 
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amount  for  transverse  displacements  on  the  order  of  a  wavelength.  The  boundary 
conditions  dictate  that  r(r)hasno  x'or  y '  dependence  at  z'  =  0,i.e.  r(r)=r(z').  Applying  the 
paraxial  approximation  to  Eq.  (13)  and  integrating  over  the  coordinate  z'  ,  one  obtains 
r(z')=a+bz'  .  The  constant  of  integration  a  is  determined  by  the  continuity  of  the  magnetic 

field  B  =  -/  (c  /  ft))  (V  X  E)  at  the  boundary.  Simplification  of  the  expression  for  yields 

\SEj.(z'  =  0)1 ,  i.e.  \SE\  as  given  by  Eq.  (12). 
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